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FOREWORD 


A  statistical  gravity  field  model  has  been 
developed  for  northern  Texas,  where  the  Gravity 
Gradiometer  Survey  System  (GGSS)  will  be  tested. 

The  model  is  described  in  this  report  using  both 
round-earth  and  flat-earth  Attenuated  White  Noise 
(AWN)  formulations.  The  model  is  fitted  to  all 
worldwide  and  local  gravity- field  information  that 
is  applicable  to  the  GGSS  test  site.  In  addition,  a 
discrete  approximation  to  the  AWN  model  is  described, 
which  can  be  used  to  compute  simulated  GGSS  survey 
data  and  surface  truth  data  for  the  test  area.  It 
is  expected  that  the  AWN  model  and  the  derived  sim¬ 
ulated  data  will  be  useful  for  optimizing,  testing, 
and  validating  GGSS  data  processing  algorithms  and 
software . 
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1. 


INTRODUCTION 


This  report  was  prepared  in  response  to  the  need  for 
a  statistical  gravity  model  that  is  tuned  to  the  north  Texas 
area  where  the  Gravity  Gradiometer  Survey  System  (GGSS)  will 
be  tested.  New  round-earth  and  flat-earth  Attenuated  White 
Noise  (AWN)  gravity  models  (Ref.  1)  have  been  developed  and  are 
described  for  the  GGSS  test  area.  These  models  are  consistent 
with  world-wide  gravity  models  and  with  surface  gravity  data 
from  the  test  site  provided  by  the  Defense  Mapping  Agency  (DMA). 
In  addition,  an  approach  is  described  for  generating  simulated 
GGSS  airborne  survey  and  surface  truth  data  for  the  Texas  test 
area.  This  approach  is  based  on  an  accurate  discrete  approxima¬ 
tion  to  the  continuous  flat-earth  AWN  gravity  model  for  northern 
Texas.  The  simulated  data,  while  different  from  actual  gravity 
values  in  the  test  area,  nevertheless  are  consistent  with  a  sta¬ 
tistical  description  of  the  earth's  true  gravity  field. 

The  discrete  approximation  to  a  flat-earth  AWN  model 
consists  of  a  set  of  mass  density  doublets,  which  are  uniformly 
spaced  on  horizontal  rectangular  grids  located  beneath  the 
earth's  surface.  (A  mass  density  doublet  is  the  mathematical 
limit  of  a  point-mass  dipole  for  which  the  separation  of  the 
positive  and  negative  masses  approaches  zero  while  the  mass 
values  approach  infinity  to  keep  the  dipole  moment  finite.) 
Each  grid  of  doublets  in  the  discrete  approximation  corresponds 
to  one  white-noise  layer  of  the  AWN  model.  Each  grid  is  at 
the  same  depth  as  the  corresponding  layer  in  the  AWN  model. 

The  spacing  of  the  doublets  is  uniform  in  each  grid  and  deter¬ 
mined  by  the  depth.  The  magnitudes  of  the  doublets  are  inde¬ 
pendent  random  variables  with  variances  that  depend  on  the 
parameters  of  the  AWN  model. 


2. 


TECHNICAL  DISCUSSION 


2.1  MASS  DENSITY  DOUBLETS 

In  this  section,  formulas  are  derived  for  the  gravity 
potential,  disturbance  vector,  and  gradient  tensor  of  a  mass 
density  doublet.  These  formulas  are  the  basis  for  computing 
realizations  of  AWN  model  gravity  fields. 

As  illustrated  in  Fig.  2.1-1,  the  gravitational  poten¬ 
tial  Tpffl(x,y,z)  at  point  (x,y,z),  caused  by  a  point  mass  singu¬ 
larity  located  at  depth  d^  beneath  the  earth's  surface  at  point 
(x,j  .y^.-dj.) ,  Sive"  by  the  formula 

* 

T^k(x,y,z)  =  ^  (2.1-1) 

Amplitude  parameter 

Distance  between  point  mass  and  observation 
point  (x,y,z) 

rijk  =  7(x-Xj)2  +  +  (z+di)2  (2.1-2) 

The  potential  T(x,y,z),  caused  by  a  mass  density  doublet 
that  is  aligned  with  the  z  axis  and  located  at  (Xj,yk,-d^),  is 
proportional  to  the  derivative  of  the  potential  of  the  point  mass 

Tljk(x,y,z)  =  c  T^k  (x,y,z)  (2.1-3) 

In  Eq.  2.1-3,  the  constant  c  =  1  km  when  z  is  measured  in  kilo¬ 


where 


meters . 
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Figure  2.1-1  Mass  Singularity  and  Cartesian  Coordinate  System 


T^k(x ,y,z)  =  -  aiJlls(Z*d-) 

rijk 


(2.1-4) 


In  Eq.  2.1-4,  the  parameter  a^j^  has  units  of  mGal'km. 

The  three  components  of  the  gravitational  disturbance 
caused  by  the  doublet  are  obtained  from  Eq .  2.1-4  by  differen¬ 
tiating  with  respect  to  x,y,  and  z: 


iik  3  aiik  (x-x-)(z+d.) 

T"Jk(x,y,z)  =  - ljk-r  — J - — 

ijk 


(2.1-5) 


3  aiik  <y-yk)(z+di> 


T?,Jk(x,y>z)  = 


(2.1-6) 


(2.1-7) 


T^Cx.y.z)  =  a^[3(fdt>2  ~  riJk] 


The  elements  of  the  gradient  tensor  (T  ,  T  ,  .  .  .  ,  T  )  are 

xx  xy  zz 

obtained  from  Eqs .  2.1-5  to  2.1-7  by  further  differentiations: 


=  .  3  auk  D^-xj)2  -  *j1k_ 


(z+di ) 


(2.1-8) 


Txyk(x’y’2)  = 


15  aiik  (x‘X.)(y-yk)(z-Hd.) 


(2.1-9) 


rijk 


T^k(*,y.*>  = 


3  «ljk  (x-x.)  [su+d^2  ■  rfik 


r1Jkl 

yy 


Tjfcx.y.z)  = 


T3Jk(x,y,z)  = 


3  aiik 

_5<y-yk>2  -  rfiR_ 

(z+dt) 

3  aijk 

< 

<y-yk> 

jk 

.5<k+di)2  -  r?iJ 

lliik 

r 

(z+di ) 

7  ' 

ijk 

[s(z+d^)2  -  3r2 -R] 

(2.1-10) 


(2.1-11) 


13) 


rijk 


2.2  AWN  GRAVITY  MODELS 


Round-earth  and  flat-earth  AWN  models  are  given  in 
Ref.  1  for  modeling  world-wide  average  gravity.  To  support 
the  GGSS  test  program,  new  round-  and  flat-earth  AWN  models 
have  been  developed  and  are  described  in  this  report.  They 
are  based  on  global  degree-variance  data  and  on  local  surface 
gravity  data  from  the  GGSS  test  area  in  northern  Texas.  The 
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!  parameters  of  these  AWN  models  are  listed  in  Tables  2.2-1  and 

2.2-2.  (In  these  tables,  D(i)  is  the  depth  of  the  i^  white- 
noise  shell  or  layer,  and  a^,(i)  is  the  rms  potential  at  the 
earth's  surface  caused  by  this  shell). 

| 

TABLE  2.2-1 


ROUND -EARTH 

AWN  MODEL 

PARAMETERS 

i 

D(i) 

oT(  i  ) 

SHELL  NO. 

DEPTH 

(km) 

rms 

POTENTIAL 
(mGal •  km ) 

1 

2.1 

2.3 

2 

5.0 

11 

3 

16 

72 

4 

52 

580 

5 

160 

2300 

6 

840 

7000 

7 

2000 

33000 

TABLE  2.2-: 

) 

FLAT -EARTH 

AWN  MODEL 

PARAMETERS 

i 

D(i) 

o^(  i  ) 

LAYER  NO. 

DEPTH 

!  (  km ) 

1 

rms 

POTENTIAL 
(mGal • km ) 

1 

l 

2  .  1 

2.3 

2 

5.0 

11 

3 

16 

72 

u 

C  O 

580 

5 

lbi 

1  2300 

i 

4 

861 

7000 

7 

2150 

33000 

-4 


The  AWN  models  in  Tables  2.2-1  and  2.2-2  were  fitted 
to  the  following  data: 

•  Gem  10B  undulation  degree  variances 
through  degree  36 

•  Rapp  180  undulation  degree  variances 
through  degree  100 

•  Along-track  surface  gravity  data  from 
northern  Texas  provided  by  the  Defense 
Mapping  Agency  (DMA). 

The  four  shallowest  white-noise  shells  of  the  new 
round-earth  AWN  model  were  fitted  to  the  northern  Texas  surface 
gravity  data  provided  by  DMA  (Ref.  2).  The  depths  and  rms 
potentials  of  these  shells  were  selected  to  yield  a  great- 
circle  gravity  disturbance  power  spectrum  that  accurately 
approximates  the  along-track  disturbance  power  spectrum  that 
was  estimated  by  TASC  using  the  DMA  surface  data.  (The  spec¬ 
trum  estimate  is  based  on  a  technique  for  state-space  modeling 
(Ref.  3),  which  allows  multiple  tracks  of  data  to  be  pooled 
for  low-frequency  spectrum  estimates  that  are  limited  by  the 
total  extent  of  the  data,  rather  than  the  lengths  of  individual 
tracks.)  The  parameters  of  the  two  deepest  shells  were  selected 
to  model  the  Gem  10B  (Ref.  4)  and  Rapp  180  (Ref.  5)  undulation 
degree-variance  data.  The  third  deepest  shell  was  selected  to 
fit  simultaneously  the  along-track  disturbance  spectrum  of  the 
northern  Texas  data  and  the  Rapp  180  undulation  degree  variances 
to  degree  100. 

Once  the  parameters  of  the  new  round-earth  AWN  model 
had  been  determined,  the  corresponding  flat-earth  AWN  model 
was  defined  as  follows.  The  flat-earth  AWN  model  has  one  layer 
for  each  shell  in  the  round-earth  model.  The  rms  undulation 
caused  by  each  layer  in  the  flat-earth  approximation  is  required 
to  be  the  same  as  the  rms  undulation  of  the  corresponding  shell 


in  the  round-earth  model.  The  depth  of  each  layer  in  the  flat- 
earth  model  is  selected  to  produce  the  same  rms  gravity  disturb¬ 
ance  at  the  earth's  surface  as  the  corresponding  shell  produces 
in  the  round-earth  model. 

The  round-earth  AWN  undulation  degree-variance  spec¬ 
trum  is  compared  with  the  corresponding  spectra  of  Gem  10B  and 
Rapp  180  in  Fig.  2.2-1.  The  AWN  model  provides  an  accurate 
fit  except  for  degrees  0,  1,  and  3.  The  AWN  model  has  signif¬ 
icant  variance  in  degrees  0  and  1,  which  is  absent  in  Gen  10B 
and  Rapp  180.  As  compensation  for  th  s  excess  variance,  the 
AWN  variance  at  degree  3  is  less  than  Gem  10B  and  Rapp  180. 

The  net  result  is  that  the  AWN  model  has  a  total  rms  undula¬ 
tion  (from  all  degrees)  of  34.5  m,  which  is  similar  to  the 
Gem  10B  rms  undulation  of  30.3  m.  Further  attention  to  details 
of  the  fit  at  degrees  0,  1,  and  3  was  not  pursued  because  the 
total  variance  of  these  terms  is  appropriate  for  modeling  the 
gravity  field  within  the  finite  GGSS  test  area  (approximately 
500  km  by  500  km  in  extent). 

The  round-earth  AWN  undulation  degree-variance  spec¬ 
trum  is  compared  with  Rapp  180  in  Fig.  2.2-2.  The  two  spectra 
are  in  general  agreement  except  for  degrees  exceeding  approxi¬ 
mately  100.  This  departure  of  the  AWN  model  from  Rapp  180  is 
the  consequence  of  requiring  the  AWN  model  to  agree  with  the 
surface  gravity  disturbance  data  from  northern  Texas. 

The  accuracy  with  which  the  flat-earth  AWN  model  approx 
imates  the  round-earth  AWN  model  is  verified  in  Fig.  2.2-3, 
which  compares  the  along-track  gravity  disturbance  spectrum  of 
the  flat-earth  model  with  the  (suitably  scaled)  great-circle 
disturbance  spectrum  of  the  round-earth  model.  As  expected, 
the  only  discernable  differences  between  the  two  spectra  occur 
at  long  wavelengths  (A.  >  2500  km). 
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Undulation  Degree  Variances  of  Gem  10B, 
Rapp  180,  and  Round-Earth  AWN  Model 
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Figure  2.2-2 


Undulation  Degree  Variances  of 
Rapp  180  and  Round- Earth  AWN  Model 
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Figure  2.2-3  Vertical  Gravity  Disturbance  Spectra  of 

Round-Earth  and  Flat-Earth  AWN  Models 


The  accuracy  with  which  the  round-  and  flat-earth  AWN 
models  approximate  the  estimated  along-track  gravity  disturbance 
spectrum  of  the  DMA  data  from  northern  Texas  is  indicated  in 
Fig.  2.2-4.  The  high-frequency  limit  of  this  comparison  is 
determined  by  the  sampling  interval  of  2.8  nm  that  was  used  in 
analyzing  the  DMA  data  sets. 

In  Fig.  2.2-4,  the  estimated  spectrum  of  the  Texas 
data  is  seen  to  flatten-off  at  frequencies  less  than  approxi- 
mately  10  cyc/km  (A  >  1000  km).  This  flattening  is  consistent 
with  the  finite  extent  of  the  local  gravity  data.  Below  the 
frequency  10  cyc/km  (wavelengths  longer  than  1000  km),  the 
AWN  spectrum  rises  so  that  the  round-earth  AWN  model  is  con¬ 
sistent  with  the  Rapp  180  undulation  degree-variance  spectrum 
for  degrees  less  than  100. 
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Figure  2.2-4 


Along-Track  Disturbance  Spectra  of  Flat-Earth 
AWN  Model  and  Texas  Data  Model 


2.3  DISCRETE  APPROXIMATION  TO  AWN  MODEL 

A  flat-earth  AWN  gravity  model  is  described  by  the 
following  parameters: 

n  The  number  of  white-noise  layers 

D(i)  The  depth  (km)  of  the  i^  layer,  for  i=l  to  n 

aT(i)  The  rms  surface  gravity  potential  (mGal'km) 
a  th 

of  the  i  layer,  for  i=l  to  n. 


The  discrete  approximation  to  the  AWN  model  is  described  using 
the  additional  parameters  defined  below: 


L  The  linear  extent  (km)  of  the  square  geographi 
area  (LxL)  over  which  the  gravity  field  is 
to  be  represented 

L(i)  The  linear  extent  (km)  of  the  i1"*1  grid  of 
doublets,  for  i=l  to  n 

s(i)  The  grid  spacing  (km)  between  adjacent  mass 

density  doublets  in  the  i^  layer,  for 
i=l  to  n 

2 

a(i)  The  rms  amplitude  (mGal'km  )  of  the  doublets 
in  the  i11*1  layer,  for  i  =  l  to  n. 


2.3.1  Parameter  Values 


Numerical  experiments  with  discrete  approximations  to 
the  flat-earth  AWN  model  lead  to  the  following  working  formulas 
for  the  doublet  grid  spacing  s(i),  the  mass  density  doublet 
rms  amplitude  ct( i),  and  the  linear  extent  L ( I )  of  the  i1*1  layer 
(for  i=l  to  n,  where  n  is  the  number  of  AWN  layers): 

s ( i )  =  0.4  D(i)  (2.3-1) 

a ( i )  =  0.3178  D(i)2  aT(i)  (2.3-2) 

L(i)  =  L  +  11.5  D(i)  for  computing  T  ,  T  ,  T 

x  y  z 

at  the  earth's  surface  (2.3-3) 

L ( i )  =  L  +  7.5[D(i)+H]  for  computing  T  ,  T  ,  ...  T 

at  height  H  xx  xy  zz 

(2.3-4) 

According  to  Eq .  2.3-1,  the  grid  spacing  is  40%  of 
the  layer  depth.  This  spacing  is  close  enough  so  that  the 
variance  of  the  gradient  field  T  changes  from  place  to  place 
within  the  LxL  region  on  the  surface  of  the  earth  by  less  than 
1%  of  its  average  value.  In  other  words,  the  40%-of-depth 
spacing  is  appropriate  to  provide  an  accurate  approximation  of 
a  stationary  gradient  field. 


The  formula  for  rms  doublet  amplitude  in  Eq.  2.3-2  is 
appropriate  for  a  grid  spacing  that  is  40%  of  the  depth.  The 
coefficient,  0.3178,  which  was  established  by  precise  numerical 
covariance  calculations,  is  close  to  the  estimated  value  of 
0.3192  that  was  obtained  from  an  approximate  analytical  co- 

t*  h 

variance  analysis.  The  linear  extent  of  the  i  doublet  grid 
specified  in  Eq.  2.3-3  is  appropriate  for  computing  vertical 
gravity  disturbance  (T  )  values  within  the  LxL  region  at  the 
earth's  surface  with  an  rms  accuracy  that  is  better  than 
0.1  mGal .  The  slightly  smaller  grid  extent  given  in  Eq.  2.3-4 
yields  gravity  disturbance  gradient  (T  )  values  with  an  rms 
accuracy  that  is  better  than  0.1  E  (0.01  mGal/km) . 

2.3.2  Computational  Formulas 

In  this  section,  formulas  are  given  for  using  an  n- layer 
AWN  model  to  compute  gravity  disturbances  and  gradients  at  an 
arbitrary  point  (x,y,z)  that  is  either  on  (z=0)  or  above  (z=H) 
the  LxL  domain.  The  calculations  consist  of  computing  n  arrays 
of  zero-mean  independent  random  variables  with  specified  vari¬ 
ances.  There  is  one  array  for  each  layer  in  the  AWN  model. 

Let  A(i)  denote  the  array  for  the  i^  layer,  and  let  A.,(i) 

th  J*' 

denote  a  number  in  its  j  row  and  k  column.  The  elements 

of  A(i)  are  computed  as  follows.  First,  a  zero-mean  unit- 

variance  random  number  v.,(i)  is  computed  for  each  grid  point 
th  J K 

in  the  i  layer.  Second,  these  numbers  are  scaled  by  the  rms 
doublet  amplitude  a(i)  (defined  by  Eq.  2.3-2)  and  stored  in 
the  corresponding  array  element: 

Ajk(i)  =  o(i)  vjk(i)  (2.3-5) 


The  expected  values  of  the  random  numbers  satisfy  the  following 
equations  for  all  admissable  i,  j,  and  k: 


(2.3-6) 


E[vjk(i))  =  0 
E[vjk(i)2)  =  1  (2.3-7) 

E[Vjk(i)vpq(m)]  =  0  for  (i,j,k)  /  (m,p,q)  (2.3-8) 

The  dimensions  of  the  n  arrays  A(l)...A(n)  are  deter¬ 
mined  by  the  linear  extent  and  grid  spacing  for  each  of  the  n 
layers.  The  number  of  elements  [ N ( i ) ]  in  array  A(i)  is  computed 
using  Eqs.  2.3-1  and  2.3-4  and  the  INT(x)  function,  which  re¬ 
turns  the  integer  part  of  x: 


N(i)  =  1NT 


L  +  7.5  [D(i )  +  H] 
- OToHl 


2 


(2.3-9) 


As  an  example  of  the  use  of  the  discrete  approximation 

to  the  AWN  model,  formulas  are  given  for  computing  the  T  (x,y,z) 

xy 

element  of  the  gradient  tensor  at  any  point  (x,y,z),  where  0  £ 
x  <  L,  0  <  y  <  L,  and  z  =  H. 


Txy(x,y,H) 


n 


2 
i  - 1 


22 

j ,keW(i,x,y) 


(2.3-10) 


In  Eq.  2.3-10,  the  limits  of  summation  on  j  and  k  are  denoted 
W(i,x,y),  which  is  square  data  window  centered  on  the  point  (x,y) 
as  depicted  in  Fig.  2.3-1.  The  sides  of  this  window  have  length 
7.5[D(i)+H]  when  computing  gradients  at  height  H.  (In  contrast, 
the  sides  of  the  window  should  have  length  11.5D(i)  when  comput¬ 
ing  gravity  disturbances  at  the  earth's  surface. )  The  window 
is  used  to  keep  computation  time  within  reasonable  bounds  by 
avoiding  a  complete  summation  over  the  entire  set  of  doublets 
for  each  observation  point.  Distant  doublets,  which  do  not 
contribute  significantly  to  the  gravity  field  at  the  observa¬ 
tion  point,  are  omitted. 
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7.5  [D(i)  +  H]  FOR  GRADIENTS 

11.5  D(i)  FOR  DISTURBANCES 


Figure  2.3-1  Data  Window  Centered  on  Point  (x,y) 

In  Eq.  2.3-10,  the  gradient  caused  by  the  jk1*1  doublet 
in  the  i*1*1  layer  is  given  by  Eq.  2.1-9,  which  is  expressed 
here  using  parameters  H  and  D(i): 


iik  15  a.  .k(x-x  .)(y-yk)[H+D( I)] 

(x  ,y  ,H)  = - ^ - J7  -  ■  - - 

rijk 


(2.3-11) 


In  Eq.  2.3-11,  the  radial  distance  is  given  as 
rijk  =  ^x-x.)2  +  (y-yk)2  +  [H+D(  i )  ] 2 


(2.3-12) 


and  the  doublet  amplitude  a^jk  is  the  appropriate  random  number 
in  the  array  A(i): 


a 


ijk 


Ajk(i) 


(2.3-13) 
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To  compute  any  of  the  other  gradient  elements  (T  ,  T  ,  .  ..) 

aA  aZ 

or  any  of  the  gravity  disturbance  components  (T  ,  T  ,  or  T  ), 

x  y  z 

the  appropriate  formula  in  Section  2.1  is  used  instead  of 
Eq.  2.3-11. 

2.3.3  Computational  Requirements 

The  amount  of  computation  that  is  required  to  compute 
simulated  GGSS  airborne  survey  data  (e.g.,  gradients)  for  the 
Texas  test  area  is  estimated  in  this  section.  The  test  area 
is  assumed  to  be  a  300-km  by  300-km  square,  and  the  airborne 
survey  is  assumed  to  have  a  height  of  0.6  km.  Therefore, 

L  =  300  km,  H  =  0.6  km,  and  Eq.  2.3-4  gives  the  following  ex¬ 
pression  for  the  linear  extent  of  the  i^  doublet  layer: 


L ( i )  =  300  km  +  7 . 5 [ D ( i )  +  0.6  km] 


(2.3-14) 


The  number  of  doublets  in  this  layer  is  N(i),  where  from  Eq.  2.3-9 


N( i )  =  INT 


300  km  +  7 . 5[ D( i )  + 
0.4  D(i) 


0 . 6  km] 


(2.3-15) 


The  number  of  doublets  in  the  data  window  W(i,x,y)  of  Eq.  2.3-10 
is  denoted  Nw(i).  For  computing  gradient  elements  aloft,  this 
number  is 


N  ( i )  =  INT 
w 


7.5  [D(i)  +  0.6  km] 
0.4  D(i) 


(2.3-16) 


Equation  2.3-16  gives  the  number  of  times  the  quantity  T^^ 

xy 

(x,y,H)  must  be  evaluated  in  Eq .  2.3-10  for  each  i  from  1  to 

n.  This  number  varies  from  N  (1)  =  625  to  N  (7)  =  361.  There 

w  w 

are  six  such  gradient  elements  (T  ,  T  ,  ...,  T  )  that  need 

xx  xy  zz 

to  be  computed,  and  the  formulas  for  them  (Eqs.  2.1-8  to  2.1-13) 


all  have  about  the  same  complexity.  That  is,  the  evaluation 
of  each  doublet  contribution  to  the  field  at  point  (x,y,H)  in¬ 
volves  approximately  10  multiplies  and  divides.  Therefore, 
the  total  number  Nfficj  of  multiplies  and  divides  to  compute  the 
six  gradient  elements  at  a  single  point  (x,y,H)  is  the  sum  of 


the  numbers  Nfficj(i)  for  each  layer: 


n 


Nmd  =  £  Nmd(i> 


(2.3-17) 


i  =  l 


Nmd(i)  =  &*Nw(iW0 


(2.3-18) 


The  total  number  of  multiplies  and  divides  for  computing 
the  six  gradient  elements  at  every  point  (x,y,H)  of  interest  in 
the  airborne  survey  is  denoted  N^..  It  is  estimated  by  multiply¬ 
ing  Nm<j  in  Eq.  2.3-17  by  the  total  number  of  points  Np  in  the 
survey  (one  point  every  1  km  along  tracks  300  km  long  spaced 
5  km  apart) : 


Np  =  301-61-2  =  3.67x10 


NT  '  NP  Nmd 


(2.3-19) 

(2.3-20) 


Using  the  AWN  model  parameters  in  Table  2.2-2,  Eqs.  2.3-16  to 
2.3-20  yield  the  following  estimate  for  the  total  number  of 
multiplies  and  divides: 

Nt  =  6.5*109  (2.3-21) 

If  the  computer  can  perform  10^  multiplies  and  divides  per 

a 

second  then  these  6.5*10  mathematical  operations  can  be  per¬ 
formed  in  1.8  hours.  The  arrays  of  random  numbers  A(l),  A(2), 
...  A(7)  range  in  size  from  382*382  for  A(l)  to  20*20  for  A(7) 
Using  single  precision,  the  total  amount  of  computer  memory 
required  to  store  these  arrays  is  approximately  730,000  bytes. 


3. 


SUMMARY 


This  memorandum  presents  a  new  statistical  gravity 
field  model  using  the  Attentuated  White  Noise  functional  form. 
The  model  describes  the  earth's  field  in  the  north  Texas  area 
where  the  GGSS  will  be  tested.  The  round-earth  and  flat-earth 
parameters  of  the  model  have  been  selected  to  match  the  sta¬ 
tistical  character  of  all  available  worldwide  and  local  gravity- 
field  information  that  is  applicable  to  the  GGSS  test  area. 

A  full  description  is  provided  of  equations  and  param¬ 
eters  that  can  be  used  to  realize  synthetic  GGSS  test  data  and 
ground  truth  data  from  the  model.  It  is  anticipated  that  the 
AWN  model  and  the  derived  synthetic  data  will  be  useful  to  the 
gradiometer  community  for  the  following  purposes: 


•  Optimization  of  GGSS  algorithms  for  down¬ 
ward  continuation  and  estimation  for  the 
residual  gravity  disturbance  vector 

•  Testing  and  validation  of  GGSS  data  proc¬ 
essing  software 

•  Analysis  of  GGSS  data-reduction  sensitivity 
to  differences  between  actual  field  statis¬ 
tics  and  practical  covariance  models 

used  for  data  processing 

•  Gaining  experience  in  the  logistics  of 
GGSS  data  processing  before  the  actual 
test  data  flow  begins. 
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